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Abstract. Stephan's Quintet (SQ) is a system consisting of at least four interacting galaxies which is well known 
for its complex dynamical and star formation history. It possesses a rich intergalactic medium (IGM), where 
hydrogen clouds, both atomic and molecular, associated with two starbursts (refered to as SQ A and B) have 
been found. In order to study the extent, origin and fate of the intergalactic molecular gas and its relation to the 
formation of stars outside galaxies and Tidal Dwarf Galaxies (TDGs), we mapped with the IRAM 30m antenna 
the carbon monoxide (CO) towards several regions of the IGM in SQ. In both SQ A and B, we detected unusually 
large amounts of molecular gas (3.1 x 10 9 Mq and 7 x 10 s Mq, respectively). In contrast, no significant CO 
detection was achieved towards HII regions south of the pair NGC 7318a/b despite their high Ha luminosities. 
The molecular gas is very extended in both SQ A and SQ B, over areas of between 15 and 25 kpc. The CO clouds 
seem to have otherwise different properties and may be of a different nature. The integrated CO line of SQ A is 
in particular much wider than in SQ B. Its CO spectrum shows emission at two velocities (6000 and 6700 kms -1 ) 
that are coincident with two HI lines. The strongest emission at 6000 kms -1 is however spatially offset from the 
HI emission and situated on a ridge south-east of the starburst region. In SQ B the CO emission coincides with 
that of tracers of star formation (Ha, 15 /im and radio continuum). The CO peak lies slightly offset from the HI 
peak towards a steep HI gradient. This is indicating that the molecular gas is forming in-situ, possibly in a region 
of compressed HI, with subsequent star formation. The star forming region at SQ B is the object in SQ that most 
resembles a TDG. 

Key words. Stars: formation - ISM: molecules - Galaxies: clusters: individual (Stephan's Quintet) - Galaxies: 
individual (NGC 7319, NGC 7318b) - Galaxies: interaction - Galaxies: ISM - intergalactic medium 



1. Introduction 

Stephan's Quintet (Hickson Compact Group 92; hereafter 
SQ) is one of the best studied examples of a Hickson 
Compact Group. It contains four interacting galaxies, 
NGC 7319, NGC 7318a, NGC 7318b, and NGC 7317. A 
fifth galaxy, NGC 7320, happens to be a foreground ob- 
ject. An optical image of the gr oup is presented in F ig. |Tj, 
(see also the excellent image in Arp fc Lorre 1976 ). One 
of its most striking properties is that the major part of the 
gas is in the intragroup medium, most likely the result of 
interactions in the past and present. 
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A plausible scenario for the dynamical history of SQ 
is presented by Moles et aL (1997). They suggest that 
a few times 10 8 yr ago the group experienced a colli- 
sion with NGC 7320c, a galaxy ~4arcmin to the east of 
NGC 73 19 but with a very s imilar recession velocity (6583 
kins" 1 , [Sulcntic ct al. 20"01 ) to the other galaxies in SQ. 
This collision removed most of the gas of NGC 7319 to- 
wards the west and east, and produced the eastern tidal 
tail which connects to NGC 7319. Presently, the group is 
experiencing another collision with the "intruder" galaxy 
NGC 7318b which strongly affects the interstellar medium 
(ISM) removed during the first collision. NGC 7318b has 
a recession velocity of 5765 km s , in contrast to the 
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Fig. 1. An archival V-band image from CFHT of Stephan's Quintet, showing the positions observed by us. NGC 7320 
is a foreground galaxy. The fourth member of the group, NGC 7320c lies about 4 arcmin east of NGC 7319. The 
velocities are taken from Sulcntic et al. (2001). The positions observed in CO are indicated by circles. The large circle 
shows the central (i.e. offset 0,0) position in each region and gives the size of the CO(1-0) beam. 



other members of the group which have velocities close to 
6600 km s -1 . Bulentic et al. (2001 ), in a multi- wavelength 



study of the group, confirm this scenario and suggest that 
the group has been visited twice by NGC 7320c. The first 
collision created the very faint tidal arm east of the inter- 
loper NGC 7320, whereas the second interaction produced 
the tidal arm which stretches from NGC 7319 eastwards. 

This violent dynamical history has induced star for- 
mation at various places outside the individual galaxies. 
ISOCAM mid-infrared and Ha observations have revealed 
a starburst region (object A in 



Xu et al. 199E, here 



Using HST observations, Gallagher et al. (2001 ) found 
115 candidate star clusters, most of them distributed 
among the tidal debris of SQ. From color-color diagrams 
they estimated their ages ranging from 2-3 Myr up to 
several Gyr. The distribution of ages sheds light on the 
star formation history in SQ: The youngest star clusters 
(with ages of less than 10 Myr) are found in SQ A and 
south of NGC 7318a/b, while somewhat older star clus- 
ters, with ages between 10 and 500 Myr are in the young 
tidal tail and around NGC 7319. This is consistent with 
the picture that the eastern tidal tail was produced in a 
previous interaction whereas the collision and star forma- 



after c alled SQ A) at the intersection of two faint optical tion a round SQ A is ongoing. Mendes de Oliveira et al 



arms s temming north from NGC 7318a/b. Several knots 
of star formation are visible in the tidal arm extending 
from NGC 7319 to th e east. At this position, identified as 
B by Xu et al. (199E ) and hereafter called SQ B, there is 
also mid-infrared and Ha emission, although much weaker 
than in SQ A. 



(2001 ) have carried out Fabry-Perot observations of the 



Ha emission in SQ A and south of NGC 7318a/b, from 
which they determined the velocity curves for several Ha 
emitting regions. SQ B was not covered by their observa- 
tions. They identified seven Tidal Dwarf Galaxy (TDG) 
candidates, selected as objects exhibiting a velocity gradi- 
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ent compatible with rotation and possessing an Lb/£hq 
ratio consistent with that of a dwarf galaxy. 

Abundant atomic hydrogen is present to the east of the 
3 central galaxies, around SQ A and south of NGC 7318a 
(jBhostak et al. 1984 |Williams ct al. 2002|, see their Fig. 



5). Williams et al. (2002) suggest that the atomic gas to- 
wards the east of SQ consists of two tidal features, each 
connected to one of the optical tidal arms (the old and 
the new one, see Fig. [l]). Molecular gas has been found at 



the position of SQ A by |Gao fc Xu (2000) with BIMA 



and by |5mith fc Struck (2001| ) with the NRAO 12m ra- 
dio telescope. Both the molecular and the atomic gas in 
this location present two velocity components, centered 
at about 6000 and 6700 km s -1 , which implies that they 
originate from different galaxies. Brainc ct al. (2001 ) in a 
study of the CO emission of TDGs detected 2.9xl0 8 M Q 
of molecular gas at SQ B. In order to follow up on these 
detections and to elucidate the extent, origin and fate of 
the CO in the intergalactic medium (IGM) in SQ, we em- 
barked on a single dish survey, the results of which we 
report in this paper. 

Details of our observations can be found in Section 2. 
In Section 3 we present the results of these observations, 
as well as a comparison between the molecular gas and 
emission at other wavelengths. We discuss our findings 
in Section 4 and summarize our conclusions in Section 
5. Following Williams et al. (2002[ ), who assumed a reces- 
sion velocity cT6^00Tm^^~aTid — ffo='75 km s _1 Mpc -1 , we 
adopt a distance of 85 Mpc, in which case 10" correspond 
to 4.1 kpc. 
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Fig. 2. CO(1-0) emission in region A, with a velocity reso- 
lution of 10.6 km s _1 . The x-scale in the individual spectra 
is velocity and ranges from 5800 to 6300 km s -1 . The y- 
scale is in T m b and ranges from -13 mK to 25 mK. The 
noise is on average 3 mK. 



2. Observations 

We observed the CO(1-0) and CO(2-l) lines at 115 and 
230 GHz in SQ in July 2001 with the IRAM 30-meter 
telescope on Pico Veleta. Dual polarization receivers were 
used at both frequencies with the 512x1 MHz filterbanks 
on the CO (1-0) line and the autocorrelator or 256 x 4 
MHz filterbanks on the CO (2-1). The observations were 
done in wobbler switching mode with a wobbler throw 
of 120" in azimuthal direction. Pointing was monitored 
on nearby quasars every 60 - 90 minutes, the rms offset 
being ss 3". At the beginning of the observations and ev- 
ery time the frequency was changed, the frequency tuning 
was checked by observing Sagittarius B2. System temper- 
atures were generally quite good, 150-200K at 115 GHz 
and 200-300 K at 230 GHz on the T£ scale. The IRAM for- 
ward efficiency, F c g, was 0.95 and 0.91 at 115 and 230 GHz 
and the beam efficiency, B c g, was 0.75 and 0.54, respec- 
tively. Our half-power beam size was 21" at 115 GHz and 
11" at 230 GHz. All CO spectra and luminosities are pre- 
sented on the main beam temperature scale (T m b) which is 
defined as T m b = {F c g/ B c g) x T£. For the data reduction, 
the spectra were summed over the individual positions and 
a constant continuum level was subtracted. 

We observed three different regions of Stephan's 
Quintet: (i) the region around starburst SQA, where 
strong mid-infrared and Ha emission had been detected 



(Xu et al. 1999), (ii) the area around a star forming re- 
gion in the eastern tidal arm of NGC 7319 (SQ B), as well 
as one position (SQ B-North) at the northern tip of the 
HI cloud covering SQ B ( |Williams ct al. 2002| ), and (iii) 
several positions to the south and north of NGC 7318a/b, 
observed by Mendes de Oliveira et al. (2001) in Ha. We 
call these positions SQ Gn where n is the number that 



Mendes dc Oliveira et al. (2001) gave to the correspond- 
ing region. The locations of our pointings are shown in 
Fig. I 

In SQ A and SQ B the spacing of the individual point- 
ings is 10" so that the CO(2-l) map is undersampled. 
Note that our observations are centered at the positions 
called SQ A and SQ B as in [Xu et a l. (1999| ) but that 



they comprise a much larger area. 



3. Results 



3.1. Molecular gas distribution 

In Fig. and Fig. |^ we present a mosaic of the CO(1-0) 
spectra observed around SQ A and SQ B. The noise in 
the spectra is about 6mK (SQ A) and 5mK (SQ B) at a 
velocity resolution of 2.66 km s _1 . The CO(2-l) emission 
was also detected at several positions where the CO (1-0) 
line was visible (see Fig. [§. The quality of the CO(2-l) 
map of SQ A is much poorer, so that no useful informa- 
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Fig. 3. CO(l-O) emission in SQ B, with a velocity resolu- 
tion of 5. 3 kms" 1 . The x-scale in the individual spectra is 
velocity and ranges from 6450 to 6750 kms -1 . The y-scale 
is in T m b and ranges from -13 mK to 38 mK. The noise is 
on average 4 mK. 



20 



-20 - 



40 



20 



-20 
offset along RA [arcsec] 



where Iqo is the velocity integrated CO line intensity ex- 
pressed in K kms -1 , D is the distance in Mpc and fl 
is the area covered by the observations in arcsec 2 (i.e. 
Q = 1.1 30 2 for a single pointing with a Gaussian beam of 
FWHM 9). We direct the reader to praine et al. (2001|) 



for a more detailed derivation of the molecular gas mass. 
When referring to molecular data from other authors, we 
scale their values to this conversion factor throughout this 
paper. 

The CO emission in SQ A is very extended. CO lines 
at a central velocity of about 6000 kms" 1 are visible in in- 
dividual spectra over an area of about 40" x 50" (Fig. ||). 
The total molecular gas mass at this velocity, derived by 
integrating over the entire area observed, is 2.2 x 10 9 M Q . 
If we were to limit the integration to the positions where 
the CO(1-0) line is clearly visible, we would obtain an 
integrated intensity of 1.2 ± 0.06 K km s -1 and a total 
molecular gas mass of 1.3 x 10 9 M Q . The error in these 
masses is about 10%, partly due to the error in the inte- 
grated intensity (5%) and partly due to the error in the 
total area which we estimate to be about 5 - 10%. Thus, 
by summing over those spectra with clearly detected sig- 
nal we apparently miss a faint, extended component which 
accounts for 40±6% of the total mass present. In the inte- 
grated CO spectrum a weaker component at 6700 kms -1 
is visible as well (Fig. ||). This component is not apparent 
in the individual spectra, and therefore we do not include 
this velocity range in Fig. |^. Its total molecular mass, in- 
tegrated over the whole area, is 8.6 x 1O 8 M . 

In SQ B the molecular gas has been detected over a 
region of about 60" x 40" (Fig. |). It is associated with 
the optically visible star forming region - the maximum 
of the CO emission coincides with it - but is more ex- 
tended. The CO emission includes the region at the very 
tip of the optical tidal tail where weaker Ha and a stellar 
condensation with blue colour (B — V = 0.2; Sulcntic et 



al. 2001) has been found. The total molecular gas mass in 
this area (7.0 x 10 s Mq , see Tab, |l|) is mor e than twice as 



much as that found by Braine et al. (2001) who observed 
only the central position. 

We did not detect CO at the majority of the positions 



of the TDG candidates observed by Mendes de Oliveira 



Fig. 4. | CO(2-l) emission in SQ B, with a velocity rcsolu- 



tion of 5.0 kms -1 . The x-scale in the individual spectra is 
velocity and ranges from 6450 to 6750 kms -1 . The y-scale 
is in T m b and ranges from -34 mK to 50 mK. The noise is 
on average 5-6 mK. 



tion can be derived from the individual spectra. Therefore 
we do not show it here. In Table |l| the integrated inten- 
sities and the derived molecular gas masses for the two 
regions are listed. All molecular gas masses in this pa- 
per are calculated using the galactic conversion factor of 
N(H 2 )/I C o = 2 x 10 20 cm " 2 (K km s" 1 )" 1 , yielding: 

M H2 [M ] = 75I co D 2 n, (1) 



et al. (2001D south and north of NGC 7318a/b. We got a 
marginal detection of Iqo = 0.5 K kms -1 when averaging 
the positions SQ C16/17 and SQ C18/19 which we used 
as an upper limit. No CO was detected in the northern 
part of the eastern HI tail (SQ B-North). This position 
was chosen because it coincides with a local peak of the 
HI emission and shows Ha emission ( [Sulentic et al. 200 1| ). 
The rms noise of the positions with non-detections are 3- 
4mK for CO(1-0) and 4-6 mK for CO(2-l), at a velocity 
resolution of 5 kms -1 . 



3.2. Kinematics and excitation 

In Fig. |5| and Fig. |6| we show the average CO(l-O) spectra 
in both regions, overlaid over the HI emission ( Williams et 
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Table 1. Molecular, atomic and ionized gas in the different regions observed 
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(1) Calculated using the Galactic conversion factor of N(H2)/Ioo = 2 x 10 20 cm -2 (K km s -1 ) -1 . 

The upper limits are calculated as a(AV/6V) 1 ^ 2 , where a is the rms noise, 8V the velocity res olution (here 5.3 km s -1 ) and 
AV the line width, adopted to be equal t o the HI line widths in this region of about 100 kms~ 1 ( |Williams et al. 2002[ ). 

(2) From HI data of [Williams et al. (2002|). 



(3) Calculated from the Ha luminosity according to SFR = 5 x 10^ a (Lu a /L e )M Q yr~ 1 ( |Hunter fc Gallagher 1986| ). 
' 4 ' Only the gas components at 6000 kms -1 . 

( 5 ) Only the gas component at 6700 kms -1 . The spatial extent cannot be quantified because the signal is too weak to show up 
in individual spectra. The value listed here refers to the area over which the line was integrated. 
' 6 ' Values for the CO and HI for the entire area with CO emission. 

' 7 ' Values for the CO and HI for the a rea within the central beam, coinciding with the Ha emission. 

^ 8 ' From Ha data of Xu et al. (1999). The first value refers to the luminosity in the starburst region, the second value to the 
total lu minosity in the region where CO was mapped 



( 9 ) From pVIendes de Oliveira et al. (2001 ), scaled down by a factor of 2.8 in order to achieve agreement with the data by 
al. (19P 



Ku et 



(io) ^he differenc e of Mh, with resp ect to the value reported in Braine et al. (2001) is due to the different adopted distances 



and the fact that Braine et al. (2001) include the helium mass in their estimate of the mass of molecular clouds. 



al. 2002), averaged over the same region. The line shapes 
in both regions are different: In SQ A the CO line is very 
broad, with a FWHM of 60 — 80 kms" 1 and line wings 
extending out to more than 200 kms -1 , showing that in 
this region large velocity gradients are present, whereas 
the line width in SQ B is much narrower, with a FWHM 
of 30 -40 kms" 1 . 

The HI gas of SQ A has two distinct velocity com- 
ponents, one at 6000 and another at 6700 kms -1 . The 
CO and HI lines generally agree well, with the exception 
of a slight velocity offset between the low-velocity com- 
ponent of the HI line, centered at 6000 kms -1 , and the 
CO line, centered at 6030-6040 kms -1 . Nevertheless, the 
line width, shapes and the relative intensity of the com- 
ponents at 6000 and at 6700 kms -1 are the same for the 
HI and CO spectra. [Smith fc Struck (200l| ) find, using 



the NRAO 12m telescope, that both velocity components 
have similar strengths. The reason for why they under- 
estimate the strength of the line at 6000 kms -1 is most 



al. 2002), and as a result is covered by their beam, but 
the 6000 kms -1 component is not. From their Fig. lg it 
is indeed obvious that their beam is smaller than the area 
mapped by us. 

The CO emission in SQ B has a central velocity of 
6625 kms -1 . There is very good agreement between the 
CO and HI lines (from Williams et al. 2002) with respect 
to the line width, central positions and special features, 
such as the blue wing. 

To better understand the dynamical behavior of the 
CO emitting gas we created several position-velocity di- 
agrams using our data. Our analysis of SQ B only gave 
marginal indications for a weak velocity gradient of about 
10 — 20 kms -1 along the east-west direction. In SQ A we 
did find evidence for a slope in velocity along a position 
angle of —40 degrees for the velocity range around 6000 
kms -1 (Fig. [?]). Moreover, for a cut perpendicular to it, 
no slope was seen. This is very similar to what we found 



likely I ecause they cover a smaller area. This can be seen 



when comparing the molecular gas masses that they derive 
for their components. Whereas their value for the com- 
ponent at 6700 kms -1 , M H2 = 7 x 10 8 M Q , agrees well 
with ours, their value for the component at 6000 kms -1 , 
Mh 2 = 7.9 x 10 s Mq, is much lower. This can be under- 
stood if the 6700 kms -1 component is more centrally con- 
centrated, as is the case for the HI emission (Williams et 



for the atomic gas using the HI data of Williams et al 
(2002|) . A similar trend has also been seen in the Ha data 



of Mcndcs de Oliveira et al. (2001) for their regions 2,3,4 
and 5. The position angle of the Ha velocity is the same 
as the one we find in CO and HI, but their observed ve- 
locity range is slightly smaller, between 5938 kms -1 and 
6010 kms -1 , whereas the CO data shows a gradient from 
about 5970 kms -1 to 6130 kms -1 . This is probably due 
to the smaller extent of only about 30" of the Ha ob- 



6 



Lisenfeld et al.: Abundant molecular gas in Stephan's Quintet 



0.005 




6000 6500 
Velocity [km/sj 

Fig. 5. CO(l-O) spectrum (full line) in SQ A, with a 
velocity resolution of 10.6 kms -1 , averaged over the total 
observed area. The dashed line represents the HI emission 
with a velocity resolution of 21.4 km s _1 in arbitrary units 
from Williams et al. (2002), averaged over the same area. 
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Fig. 6. CO(1-0) spectrum (full line) in SQ B, averaged 
over the positions with offsets between 30" and -10" in 
RA and -10" and 10" in Dec. The velocity resolution is 5.3 
km s _1 . The dashed line represents the HI emission with a 
velocity resolution of 21.4 kins -1 in arbitrary units from 
Williams et al. (2002), averaged over the same area. 

servations. Normalized to the same spatial extent the Ha 
gradient is comparable to the gradient in CO within the 
uncertainties. 

In Fig. H and || we show the CO(2-l) lines, aver- 
aged over the same area as the corresponding CO (1-0) 
lines. In order to compute the CO(2-1)/CO(1-0) inte- 



6100 



6000 - 




< 



-40 
SE 



-20 
angular offset [arcsec] 



20 



NW 



Fig. 7. Position- velocity diagram of the CO(1-0) emission 
in SQ A at a position angle of —40 degrees. The contours 
are at 4.5, 7 and 9.5 mK (T m b). 
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Fig. 8. CO(2-l) emission (full line) in SQ A, averaged 
over the total observed area. The velocity resolution is 
10.6 kms -1 . The dashed line represents the HI emission 
with a velocity resolution of 21.4 km s _1 in arbitrary units 
from Williams et al. (2002) , averaged over the same area. 



grated intensity ratio in SQ A and SQ B, the CO(2-l) 
spectra have been convolved to the spatial resolution of 
the CO(1-0) beam. To minimize the error of this ratio, we 
have included only spectra at positions with a S/N > 3er 
of the CO(l-O) line. In SQ A, the CO(2-1)/CO(1-0) ratio 
is 0.69 ± 0.16 and in SQ B 0.56 ± 0.13. Within the errors, 
there is no significant difference in the line ratios. They are 
slightly lower than the average found for a mixed sample of 
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Fig. 9. CO(2-l) emission (full line) in SQ B, averaged 
over the same area as in Fig. [| and with a velocity res- 
olution is 5.0 kms" 1 . The dashed line represents the HI 
emission with a velocity resolution of 21.4 kms -1 in ar- 
bitrary units from Williams et al. (2002), averaged over 
the same area. 



nearby galaxies ( CO(2-1)/CO(1-0)= 0.89 ± 0.06, |Brainc 
& Co mbes 1992 ). Better CO (2-1) data are necessary in 
order to judge whether this difference is significant. 

3.3. Comparison of the CO and HI distribution 

In spite of the practically identical kinematics, the atomic 
and molecular gas distributions show differences both in 
SQ A and in SQ B. In SQ A, the CO lines are wide and 
have low maxima. As a result, it is difficult to display 
the velocity integrated intensities of the CO emission be- 
cause of their large uncertainties. We therefore show in 
Fig. O a comparison of the HI emission around 6000 
kms with the observed positions in CO. The positions 
where we clearly (> 3cr) detected a line are indicated by 
a cross and the observed positions as circles. Note that, 
as mentioned earlier, also at the positions without a clear 
CO detection some faint emission must be present. From 
Fig. [l0| it is clear that the overall distribution of the CO 
emission is different from the HI emission. Even though 
some CO marginally coincides with the HI peak, most of 
it is located to the south-east along a faint HI ridge. 

In Fig. |ll| we display the contours of the integrated 
CO intensity of SQ B overlayed onto contours of an HI in- 
tegrated surface brightness map. The maxima of the CO 
and HI emission are close, but do not exactly coincide. 
Instead, the CO peaks slightly (10 - 15 ") west from the 
HI ridge, towards a steep HI gradient, suggestive of com- 
pression in the gas which implies that this over-density 
in HI may have provoked the formation of molecular gas. 
The molecular-to-atomic gas ratio is highest at the peak 



P 



(from Xu et al. 
around 6000 kms" 



1999| ), covering 
, overlaid with 
levels at 



Fig. 10. A grey scale image of the Ha emission in 
SQ A obtained with a filt er centered at 66 67 A and 
a FWHM of 66 A 
the Ha emission 

an HI contour map ( Williams et al. 2002 
1.2, 2.4, 3.6, 4.8, 6.0, 7.2, 8.4 x 10 20 atoms) showing the ve- 
locity integrated column density of the line at 6000 km s~ 1 . 
The circles show the location of our pointings and the 
crosses indicate the positions where CO emission was de- 
tected above a level of 3o\ 



of the CO emission and west of it. This means, under the 
assumption that the molecular ga s is being formed i n situ, 
as has been suggested for TDGs (Braine et al. 2001), that 
this compressed region is the place of the most efficient 
transformation of HI into H2. 

The HI mass in SQ A, integrated over the area mapped 



in CO, is 1.6 x 10 9 M (9 x 10 8 M ) for the velocity com- 
ponent at 6000 kms" 1 (6700 kms" 1 ) (see Tab. |). This 
means that the molecular gas mass is comparable to the 
atomic gas mass with a ratio of Mh 2 /Mhi= 1.4 (1.0). The 
HI mass of the entire area around SQ A is 3.1 x 10 9 M , 
which is of the same order as the total mass of the gas in 
molecular form. The atomic gas mass in SQ B, integrated 
over the region where CO was detected, is 1.4 x 1O 9 M . 
This yields a molecular-to-atomic gas mass ratio of 0.5, 
higher than t he values of less than about 0.2 found for 
other TDGs ( |Braine et al. 200l| ). 



3.4. Comparison of the CO and Ha distribution 

The Ha emission in SQ A consists, as do the HI and 
CO data, of two velocity components, centered around 
6000 kms -1 and 6700 kms -1 , respectively (Sulentic et al. 
|2001| , |Iglesias-Paramo fc Vilchez 2001| , [Xu et al. 1999| ). 
The distribution of both velocity components is similar 
in the regi ons that our CO observations covered. Using 
the data of |Xu et al. (1999 ), taken with narrowband fil- 



ters, we present in Fig. |lO|the spatial distribution of the 
Ha emission at 6000 kms -1 and our CO detections. We 
observe that the CO gas follows closely the Ha emission 
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Fig. 11. A grey scale image of the Ha emission in SQ B 
obtained with a filter centered at 6737 A and a FWHM 



of 76 A (from Xu et al. 1999 ), covering the Ha emis- 
sion around 6600 kms -1 . overlaid with an HI contour 



map ( Williams et al. 2002 , dashed contours, levels at 



1.2,2.4,3.6,4.8,6.0,7.2,8.4 x 10 2 " atoms cm" 2 ) showing 
the velocity integrated column density of the line at 6600 
kms -1 , and thick (red) contours representing Ico, the ve- 
locity integrated CO emission, (contour levels at 0.4, 0.65, 
0.9 and 1.15 K km s" 1 ). 



in the north-south direction, but surprisingly little CO is 
seen close to the brightest Ha peak in the starburst region 
at SQ A. 

Similarly, we present in Fig. [ll]the Ha image of SQ B 
overlayed with the contours of the molecular gas. The Ha 
emission is situated 5" north of the maximum of the CO 
emission and approximately 15" north-west from the max- 
imum of the HI emission. We conclude that the Ha and 
CO emission coincide within the limited spatial resolution 
(21") of our data. Also, the velocities agree, with the veloc- 
ity of the CO line bein g 6625 kms' 1 and th at of the Ha 
emission 6617 kms -1 ( Sulcntic et al. 2001 ). Altogether, 
this is pointing to a causal connection between molecular 
gas and recent star formation. 

In Tab. we list the molecular gas mass divided by 
the SFR (derived from the Ha luminosity) which gives 
the gas consumption time (i.e. the inverse of the star for- 
mation efficiency). The values we derive for SQ A and 
SQ B are typical to those of spiral galaxies, which have 



been estimated to lie between 4 and 40 x 10 yr (Kennicutt 
1998()7] 

3.5. Molecular gas fraction and star formation in the 
non-detections 

We observed various positions (see Fig. [j]) south of 
NGC 7318a/b and north of SQ B in CO without getting 
a detection. The results are summarized in Tab. [|. 



The position north of SQ B (SQ B-North) coincides 
with the north tip of the HI cloud and was chosen because 
|Sulcntic et al. (2001 ) detected Ha emission at the same 



velocity as the HI and suggested that an extragalactic HII 
region may be forming. Gallagher et al. (2001) identified 
two starclusters (30 and 31) at the same position. The 
clusters have blue colours and are therefore young, less 
than 10 - 500 Myr, and must have formed in situ. At this 
position we derived an upper limit of Mh 2 /Mhi< 0.08. 
This value is much lower than at the other places in SQ 
where we detected molecular gas, bu t it is not unusua l 
compared to typical values for TDGs ( Braine et al. 2001 ). 
The upper limit of M H2 /SFR< 4 x 10 s yr is at the low 
end of the typical range for spiral galaxies. 

Around NGC 7318a/b we observed several of the po- 
sitions where Mendes de Oliveira et al. (2001) have found 
Ha emission (their positions 8, 16/17, 18/19, 20, 21, 22 
and 23; note that positions 16 and 17 as well as 18 and 
19 fall within the IRAM beam). The authors suggest that 
some of them (8, 20, 21, 22 and 23) are TDG candidates, 
mainly based on the detection of a velocity gradient in the 
Ha data which they interpret as rotation. We observed in 
addition to these TDG candidates positions 16/17 and 
18/19 because of their relatively high HI column density. 
The upper limits of the molecular gas fraction, Mh 2 /Mhi, 
are low for some objects, but not unusual compared to 
spiral galaxies and TDGs. The upper limits for the gas 
consumption time, Mh 2 /SFR, on the other hand, are, for 
mos t objects, much low er than the values found for TDGs 
(see Braine et al. 2001 ) or spiral galaxies. This is particu- 
larly surprising, as the large size of the CO(1-0) beam of 
21", corresponding to 8.6 kpc, should comprise the molec- 
ular gas associated with the star formation even if the Ha 
and CO emission do not coincide exactly, because, e.g. the 
ongoing star formation has destroyed the molecular gas lo- 
cally. A low metallicity could explain the non-detections 
but seems unlikely as we did detect CO at other places 
in SQ. We checked for the presence of star clusters from 
|Gallagher et al. (2001 ) at the observed positions in or- 



der to find out whether star clusters have been formed 
since the beginning of the collision with NGC 7318b of 
less than 10 8 yr ago (Sulcntic ct al. 2001). Region C23 



coincides with an old star-cluster (number 117), with an 
age of several Gyr, which was therefore born before the 
start of the latest interaction. Regions C8, C16/17 and 
C18/19 coincide with star clusters 133, 134, 135, 137, 122, 
128, and 130, some of which are younger than 10 Myr. At 
the other observed positions no star clusters were found 
within a radius of 10". Thus, young star clusters formed 
after the beginning of the collision have only been found 
close to the parent galaxy NGC 7318a and not in the TDG 
candidates C20 - C23. 



3.6. Comparison of the CO distribution to other tracers 

In Fig. [l2|, P and P we compare the CO distribution to 
bhat of the 15 /im emission observed by ISO CAM (Ku ct 
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Fig. 12. A grey scale optical image with a contour o verlay 
of the 15 ami emission by ISOCAM flXu et al. 1999|). The 
thick ( red) contours on SQ B and the circles on SQ A 
show our CO observations, as described in Fig. |ld| and |li|T 



Fig. 14. A grey scale optical image overlaid with contours 



of th e 20cm radio continuum emission ( Williams et al 



2002 ). The thick (red) contours on SQ B and the circles 
on SQ A show our CO observations, as described in Fig. 
Hand ID] 



from the CO maximum and is located towards the peak 
of the HI surface density. 

In t he intergalactic medi um around SQ A, ra dio con 
tinuum ( Williams et al. 2002 ) and X-ray emission ( Pietsch 



et al. 1997|) are present. As previously noted (e.g. pulcntk 



et al. 2001), there is good agreement between the radio 



J if 




Fig. 13. A grey scale optical image overlaid wit h contours 
of the X-ray HRI emission from ROSAT (see Pietsch et 
al. 199^), smoothed to a resolution of 10". The circles on 
SQ A show our CO observations, as described in Fig. |To|. 



At SQ |B, no X-ray emission was found. 



al. 1999), the X-ray emission from the ROSAT archive 



and the 20cm radio continuum emission ( Williams et al 
2002[)7|respectively. 

The overall distribution of the dust traced by the 
ISOCAM 15 jum image is similar to the CO emission 
in both SQ A and SQ B (Fig. |l2|). In SQ A the 15 fxm 
emission coincides with the maximum in the Ha and indi- 
cates the position of the intergalactic starburst. The CO 
is slightly offset towards the south and west of the 15 /im 
peak. The CO emission coincides with faint 15 /im emis- 
sion at all our observed positions, showing that the pres- 
ence of molecular gas is accompanied by the presence of 
dust. In SQ B the 15 /im peak is marginally offset (5-10") 



continuum, the X-ray and Ha emission in the elongated 
north-south feature north of NGC 7318a/b. All along this 
ridge we find CO emission (Fig. [l^ and |l4|). CO is also 
present in a region without X-ray emission towards the 
north which coincides with the edge of the ISOCAM peak. 
In SQ B no X-ray emission is detected, but weak (2a) ra- 
dio continuum emission can be seen, coinciding with the 
peak of the CO emission. 

The radio continuum emission in both regions is mostly 
non-thermal (i.e. synchrotron emission) as can be deduced 
by calculating the expected thermal radio emission from 
the Ha emission. The thermal radio emission derived from 
the H a flux yields (following eq. (3) and (4) in |Condon 
1992), a few percent of the radio continuum flux in SQ A 



and about 10% in SQ B. This thermal fraction at 1.49 
GHz is typical for star- forming galaxies (Condon 1992). 



4. Discussion 

4.1. Molecular gas, shocks and star formation in SQ A 

The enormous amount of molecular gas that we have 
found and its large extent are surprising. In the follow- 
ing we will discuss how this molecular gas fits into the 
current picture of the IGM in SQ A. 

Th e IGM around SQ A shows two velocity compo- 
nents. Sulentic et al. (2001 ) suggest that the gas at 6700 



km s was stripped during the earlier visit of NGC 7320c 
and part of it has become shock-ionized due to the present 
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interaction with NGC 7318b. In this view, the Ha emis- 
sion at 6700 kms" 1 traces the shock front. Strong support 
for this shock hypothesis is the fact that X-ray (tracing 
hot gas) and radio continuum emission (tracing cosmic ray 
electrons produced in shocks) are coincident with the Ha 
emission. The gas at velocities around 6000 kms -1 orig- 
inates from NGC 7318b and is being removed from this 
galaxy in the ongoing collision. The Ha emissions at 6000 
and 6700 kms -1 overlap to a large extent in the region 
where we observed CO. In HI, the two velocity compo- 
nents have their maxima at the same position, but the 



distributed and therefore not detected by interferometric 
observations. 

The origin of the molecular gas in SQ A is unclear. 
It could have been stripped together with the atomic gas 
from NGC 7318b or it could have formed in situ. In the 
center of NGC 7318b, a molecular gas mass of 1.2 x 10 9 
M Q has been found ( |Smith fe Struck 200 1| ). This yields, 
together with its blue luminosity of 4.6 x 10 10 M Q , a ratio 
of L-q/Mh 2 = 38, considerably higher than the average 
value of Lb/Mh 2 as 6 (with a scatter of a factor 3) th at has 



compotient at t>YUU kms is spatially much les s e xtended ' 



than t hat a t bUUU kms' 1 ( [Williams ct al. ^UU^| ). pulcntic 
et al. ( ^001 ) suggest that the apparent spatial coincidence 
is a projection and that NGC 7318b, and thus the shock 
caused by it, are moving along the line of sight. 

The presence of molecular gas indicates a rather 
cool environment. It is therefore surprising to find large 
amounts and extended CO emission at positions where 
shocks are present and are heating and ionizing the gas. 
J-shocks with velocities above about 50 kms -1 are able 



been found for isolated galaxies of this luminosity ( Perea 
ct al. 1997| ). Thus, the molecular gas in SQ A could have 



to effectively destroy molecular gas (Hollenbach & McKee 
1980|~ f989| ). Shocks with these velocities are plausibly 



present given the large velocity difference (600 - 1000 
kms -1 ) between the intruder NGC 7318b and the rest 
of the SQ galaxies. 

This strongly suggests that shocks do not affect the 
gas at 6000 kms -1 , the velocity at which we found the 
major part of the CO emission. This means that the X- 
ray and radio continuum emission tracing shocks are only 
associated with the gas at 6700 kms -1 . In addition to 
the presence of extended molecular gas, support for this 
hypothesis comes from the atomic gas. Similarly to the 
CO, the HI component at 6000 kms -1 is very extended, 
indicating that shock heating is not affecting this area so 
that cool gas can exist. Contrary to this, the HI compo- 
nent at 6700 kms -1 is spatially much more concentrated. 
Furthermore, the spatial coincidence between the 15 /mi 
and and CO emission support the absence of shocks be- 



cause dust is very easily destroyed by them (e.g. Borkowski 
& Dwck 1995, Jones et al 



199§ ) 



The molecular gas mass in the starburst region^] is 
7.0 ± 1.0 x 10 8 Mq. The surprising result of our obser- 
vations is that the molecular gas content of the starburst 
region is only a small part (about one fourth) of the to- 
tal molecular gas content and that the starburst does not 
mark the maximum of the CO emission. Our molecular gas 
mass in the starburst region is only slightly higher than 
the mass estimated from interferometric observations by 
Gao fc Xu (2000D (2.4 - 5.3 x 10 s M ). The observa- 



tions of Gao & Xu (2000) showed CO emission only at 



the position of the starburst. This could indicate that the 
extended emission that we have found is more smoothly 



1 I.e. adding the 4 pointings around the position (-5", 5"), 
where the maximum of the Ha emission is found, and summing 
both velocity components. 



its origin in NGC 7318b. On the other hand, if the molec- 
ular gas has formed in situ, a large scale process and not 
gravitational collapse must be responsible for its forma- 
tion over a scale of about 25 kpc. The velocity of this gas, 
6000 kms -1 (which is higher than the central velocity of 
NGC 7318b of 5765 kms -1 ) could be an indication that 
it has already been accelerated by the collision with the 
intergalactic gas and the slight offset between the molec- 
ular (6030 kms -1 ) and atomic gas (6000 kms -1 ) central 
velocities might show that the molecular gas is forming in 
the most accelerated parts. 



4.2. Origin of the molecular gas in SQ B 

The molecular gas in SQ B could have either been stripped 
from the parent galaxy together with the atomic gas or it 
could have formed in situ. Studying a sample of TDGs, 
Braine ct al. (200l| ) argue in favor of the in-situ forma- 
tion. Their main argument is that the CO distribution in 
spiral galaxies is normally concentrated towards their cen- 
tral parts whereas the atomic gas is much more widely dis- 
tributed and therefore more easily stripped. In SQ the sit- 
uation is more complicated because more than two galax- 
ies are interacting. However, the arguments of Braine et 
|al. (2001[ ) seem to be applicable to NGC 7319, the galaxy 
thought to be the progenitor of the eastern HI cloud and 
the young optical tidal tail. All of the HI from this galaxy 
has been stripped, whereas CO is still found within its 
optical disk ( |Yun et al. 1997| , |Smith fc Struck 200l|). Its 



molecular gas content of 4.8 x 10 M Q (Smith & Struck 



200 1|) and its blue luminosity (Lb = 5.5 x 10 10 L ; |Verdcs- 



Montcncgro et al. 1998| ) yield Lb /Mh 2 = 11 which is nor 



mal compared to other galaxies (see above in Sect. 4.1). 
This supports the idea that in NGC 7319, as in other inter- 
acting galaxies, the HI was more easily stripped whereas 
CO stayed within the optical disk. 

In SQ B the molecular gas is found in close association 
with the atomic gas and with indicators of star formation, 
as described in detail in the previous sections. Altogether, 
this spatial and (in the case of HI also kinematical) co- 
incidence suggests that we are seeing in-situ molecular 
gas formation possibly provoked by compression of the 
atomic gas, and the subsequent star formation traced by 
Ha emission and other tracers. Further support for this 
comes from the distribution of Mh 2 /Mhi- If the molecu- 
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lar gas had been stripped together with the atomic gas, 
the molecular-to-atomic gas mass ratio would have been 
preserved, il no molecular gas formation or destruction 
took place after the stripping. In this case it would be 
an unlikely coincidence that we find the highest value of 
Mh 2 /Mhi close to the HI peak. 



4.3. Tidal Dwarf Galaxies in Stephan's Quintet? 

We observed very different regions in SQ and our CO 
observations have shown that the properties of the ISM 
are very diverse. Some of our observed regions have been 



claimed to be TDG candidates (e.g. 


Mendes dc Oliveira 


et al. 2001 




Iglesias-Paramo & Vi'lchez 2001 


)■ 


Due et al. 



(2000) defined a TDG as a self-gravitating entity formed 
out of tidal debris. Although it is, without a detailed kine- 
matical study, difficult to decide w hether a given TDG 
candidate is indeed self-gravitating, Braine et al. (2000 ) 
have argued that the detection of molecular gas form- 
ing in situ gives strong support to the hypothesis of self- 
gravitation. Most objects observed in SQ are, however, 
different in one way or another f rom t he sample of TDGs 
observed by [Braine et al. (2000| ; |200l| ). 



TDG, i.e. an object that is bound by gravitation. These 
doubts are reinforced by the large dynamical mass that 
is deduced from the CO velocity gradient in SQ A (see 
Sect. 3.2). Under the assumption that the region around 
SQ A is gravitationally bound and that the velocity gradi- 
ent is due to rotation we can calculate the dynamical mass 
(Mdyn = RAV 2 /G, where R is the radius of the region, 
AV is the velocity gradient along R and G is the gravita- 
tional constant). With R = 30" and AV = 80 kms" 1 we 
derive M dyn = 1.9 x 10 10 M Q . This is much larger than 
the mass of the molecular and atomic gas at this velocity 
(3.8 x 10 9 M ). The st ellar mass in the starburst region 
has been estimated by Xu et al. (1999 ) from a K-band 
image to be 0.8 — 1.6 x 10 7 M©, much smaller than the 
gas mass. Thus, the dynamical mass is much higher than 
the visible mass meaning that this region would need a 
high fraction of dark matter in order to be gravitationally 
bound. This is in contrast to what has been found for a 



SQ B is the object that most closely resembles TDGs. 
Its has a clear tidal origin and we have summarized evi- 
dence that the molecular gas is forming in situ. The ratio 
of molecular gas mass to the Ha luminosity is normal for 
star forming galaxies. The star formation rates that can 
be derived from different tracers such as Ha emission, 15 
/zm emission (following the estimate in 



Xu et al. 



1999| ) 

or radi|o continuum (applying eq. 'li m Condon (lUij^ ), 
extrapolated with a Salpeter initial mass function down 
to stars with masses above 0.1 M Q ) give similar values of 
about 0.2 M Q yr _1 . These values are in the same range 
as those found for TDGs (0.01 — 0.3 Mayr" 1 ) and severa l 
nearb y dwarf galaxies (0.05 — 0.9 M Q yr _1 ) ( Braine et al 



sample of TDGs ( Braine et al. 2001 ) where the visible and 
dynamical masses were found to be very close obviating 
the need for dark matter. We interpret the high dynami- 
cal mass as evidence that this region is not gravitationally 
bound on the size scales observable at the current resolu- 
tion. The extended and homogeneous distribution of the 
molecular gas in this region makes it unlikely that a gravi- 
tational collapse is responsible for its formation. The star 
formation in the starburst region SQ A has rather been 
triggered by local processes, as e.g. by the pressure of the 
hot surrounding gas, as suggested by (Xu et al. 199E). 

Most of the objects in SQ C possess a gradient in the 
Ha emission indicative of rotation (Mendes de Oliveira 
et al. 2001] ) and might therefore be self-gravitating ob- 



200l| ) showing that the observations of SQ B are consis- 
tent with the formation of a dwarf galaxy in the tidal tail. 
However, the large extent of the molecular gas covering 
the entire tip of the optical tail is surprising and different 
from other TDGs. One reason is probably the fact that 
star formation happens at two positions along the optical 
tidal tail: apart from SQ B there is a weaker star-forming 
knot showing a blue colour and Ha emission at the very 
end of the optical tidal tail (at 02000 = 22 h 36 m 12.2 s and 
^2000 = 33°57'30"). Molecular gas was detected at this 
position. Higher resolution observations are necessary in 
order to find out whether the molecular gas is concen- 
trated towards the two star forming regions as would be 
expected if TDGs are forming. 

The situation in SQ A is entirely different. Whereas 
the star formation is very concentrated at the central po- 
sition of SQ A, the distribution of the molecular gas is 
much more extended with no maximum at the starburst. 
This lack of association of the molecular gas with the star- 
burst region, together with the unusual origin of the ISM 
in this region brings up strong doubts whether SQ A is a 



jects. However, the lack of gas, both atomic (except for 
C16/17 and C18/19 the atomic gas densities are < 10 20 
cm~ 2 ) as well as molecular, makes them different from 
other TDGs. For this reason it is unlikely that they will 
become like TDGs that resemble dwarf irregulars or blue 
compact dwarfs as they do not have the gas reservoir to 
continue forming stars for a long time. 



5. Conclusions 

We have observed CO in several regions in SQ. Our main 
results are: 

1. In both SQ A and SQ B we detected large amounts 
of molecular gas of 3.1 x 10 9 M and 7.0 x 10 s M , 
respectively. The molecular-to-atomic gas mass ratios 
are high (1.2, respectively 0.5), much larger than what 
has been found in typical TDGs (Braine et al. 2001). 
The molecular gas is very extended, over regions of 
about 15 - 25 kpc. 

2. The CO spectrum in SQ A consists of lines at differ- 
ent velocities, centered at 6030 and 6700 kms -1 , with 
most of the emission coming from the low velocity com- 
ponent, whereas in SQ B the CO is at 6625 kms -1 . 
The kinematical agreement (velocity, width and spec- 
tral shape) between the CO and HI lines is very good 
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in both regions, except for a small offset in the low- 
velocity components of the CO and HI lines in SQ A 
(the HI line peaks at 6000 kms" 1 ). 
In SQ A, the CO emission at 6000 kms -1 is found 
south of the starburst region, covering the southern 
part of the HI distribution at this velocity and coin- 
ciding with X-ray and radio continuum emission which 
are indicative of shocks. The much weaker CO emission 
at 6700 kms -1 is more spatially concentrated. We con- 
clude that the presence of large amounts of extended 
molecular gas at 6000 kms -1 implies that the shocks 
traced by X-ray and radio continuum emission only 
affect the gas at 6700 kms -1 and that the apparent 
spatial coincidence is a projection. 
In SQ B there is good spatial agreement between the 
tracers of star formation (Ha, 15 /im and radio contin- 
uum) and the CO emission. The CO peak is slightly 
offset from the HI peak towards a steep HI gradient. 
We interpret this as indicating that the molecular gas 
is forming in-situ, possibly in a region of compressed 
HI and that the star formation is taking place subse- 
quently. This, together with the fact that the emission 
at different wavelengths is similar to that of normal 
star forming galaxies and that SQ B has a clear tidal 
origin, makes it the object in SQ that most resembles 
a TDG. 

We have searched without success for CO in several 



objects around NGC 7318a/b in which Mendes de 



Oliveira et al. (2001) found Ha emission. The up- 



per limits for the ratios between molecular gas mass 
and SFR (derived from the Ha emission) are lower 
than those found in spiral galaxies and TDGs. The 
molecular-to-atomic gas mass ratio is also low, though 
not unusual. The low atomic and molecular gas con- 
tent of these regions makes them different from TDGs 
and will not allow them to sustain star formation for 
a long time. 
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